Tightly Localized Stationary Pulses in Multi-Level Atomic System 
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We show the pulse matching phenomenon can be obtained in the general multi-level system with 
electromagnetically induced transparency (EIT). For this we find a novel way to create tightly 
localized stationary pulses by using counter-propagating pump fields. The present process is a 
spatial compression of excitation so that it allows us to shape and further intensify the localized 
stationary pulses, without using standing waves of pump fields or spatially modulated pump fields. 
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Recently, an important progress in electromagnetically 
induced transparency (EIT) [H, 0, i, H i, H, Q is that 
experimental realization for the coherent control of sta- 
tionary pulses was achieved by using standing waves of 
pump fields in the three- level system [§] . The creation 
of stationary pulses can enhance the nonlinear couplings 
between photons or collective excitations corresponding 
to stored photons, both of which are useful for deter- 
ministic logic operations. The key point in the creation 
of stationary pulses is expressed by the pulse matching 
phenomenon between the forward (FW) and backward 
(BW) propagating probe fields [1, @, For a three- 
level system, the technique to generate tightly localized 
stationary pulses involves the use of standing waves of 
pump fields with a frequency-comb or spatially modu- 
lated pump fields However, such tight localization 
is created by a filtering process rather than a compres- 
sion of excitation. Thus the three-level technique cannot 
be applied directly to applications in quantum nonlinear 
optics. 

On the other hand, coherent manipulation of probe 
lights has been studied in the four-level double A-type 
system ll|, |l2| and also in the general multi- level atomic 
system that interacts with many probe and pump fields 
I3L 1^ . It has also been shown in Ref. 13 1 that, one can 
convert different probe lights into each other by manip- 
ulating the external pump fields based on such general 
EIT method, indicating a sort of pulse matching phe- 
nomenon between all applied probe fields. This observa- 
tion also motivates us to probe into a new technique of 
creating localized stationary pulses based on multi-level 
atomic system. In this rapid communication, we shall 
show the tightly localized stationary pulses can be ob- 
tained through a spatial compression of excitation the 
general multi-level EIT system. 

We consider the quasi-one dimensional system shown 




FIG. 1: (color online) (a) General m-level atomic system cou- 
pled to m — 2 quantized probe and classical pump fields which 
propagate in +z or —z directions, (b) No. 1 to No. m — 3 
pump/probe pulses propagate in the -\-z direction, while No. 
m — 2 pump/probe pulse propagates in the —z direction. 



in Fig. 1(a) for an ensemble of m-level atoms interact- 
ing with m — 2 quantized probe fields which couple the 
transitions from the ground state \b) to excited state \ea) 
(1 < cr < TO — 2) with coupling constants g^, and m — 2 
classical pump fields which couple the transitions from 
the state |c) to excited ones \ea) with Rabi- frequencies 
57o.(z,t). All probe and pump fields are co- propagating 
in the -l-z or — z direction (Fig. 1(b)), and 

EAz,t) = y|^4(z,t)e^('=— *), 



(1) 



where a — 1,2,. 



2, Ea and fio-o are slowly- varying 
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amplitudes, kp„ and fcco-, respectively z-component wave 
vectors of probe and pump fields, can be positive or neg- 
ative. For fcpcr > and kca > {kpa < and fcccr < 0), it 
means the trth pair of probe and pump fields propagate in 



2 



the +z (— z) direction. We consider all transitions to be 
at resonance. Under the rotating-wave approximation, 
the interaction Hamiltonian can be written as: 

r dz 

V = - / —{hN^g,aeAz,t)£,{z,t) + 

(T = l 

m-2 

+ hNY,^<rQ{t)ae^c{z,t)+h.c), (2) 

(T=l 

where N is the total atom nmubcr, L is the length of the 
medium in the z direction, and the continuous atomic 
variables (Tfj,^{z,t) — -^^z-eN^ '^jj.ui.i) ^^'^ defined by 
a collection of iV^ 1 atoms in a very small length 
length interval Az [3]. a^^^ = jgi ) (feJ | e-'C^p"^-'^--''*) 
and ai^^ = |e^)(c'' | Q-i-ii^o„z-uia^at) g^j.g ^j^g slowly- varying 
parts of the jth atomic flip operator. We note that an es- 
sential difference between our model and the three-level 
case is that for the case of multi- frequency optical pulses, 
here the one- and two-photon detunings can be avoided 
for all optical transitions, and no standing waves of pump 
fields or spatially modulated pump fields are used. 

The evolution of the slowly- varying amplitudes £a{z, t) 
can be described by the propagation equations 



d_ 

dt 



£a{z,t) = igaNabe^iz,t), 



(3) 



where we note v^/kpa — ±c for the ±z directional prop- 
agation field. Under the condition of low excitation, 
i. e. (Tbb « 1, the atomic evolution governed by the 
Heisenberg-Langevin equations can be obtained by 



CTbe^ = —Jbe^<Jbe^ + IQaSa + i^aO<^bc + Fbe 



m-2 



m-2 



F, 



be : 



(4) 



(5) 



a=l 



For our purpose we shall calculate abc to the first 
order, neglecting the small time derivatives of f2o-0j thus 



CTbc 



1 . 1 

cr = l 
m-2 



r 



Q4 ^ 939kgMao£j£k£a 
3ka 



(7) 



The second term in the right hand side of above equa- 
tion represents the nonlinear couplings between the probe 
pulses. 

The dark-state polaritons (DSPs) in the general multi- 
level EIT system is first obtained in [3] , where the single- 
mode probe pulses are considered. Accordingly, the dark- 
and bright-state polaritons (BSPs) in the present general 
multi-level system can be defined by: 

i) = COS 6* cos(f>j£i 

m — 2 m — 3 

+ COS0 ^ sin Y[ COS (j>j£i 

1=2 j=l 

(8) 



-sin6'(<) VNabciz,t) 



m— 3 

4(z, i) = sin6' J]^ cos0j£i 

m — 2 m— 3 

-1-sinf? sm(f)i-i Y\ cos(j)j£i 

1=2 3=1 

+ cose{t) y/Nabc{z,t), 



(9) 



which are superpositions of the atomic coherence and the 
m — 2 probe fields. The mixing angles 9 and in the 
new quantum fields are defined through 



(6) 



where 7^1^ are the transversal decay rates that will 
be assumed jbe„ = F in the following derivation and 
F^i, are (5-correlated Langevin noise operators. From 
the Eq. (|4]) we find in the lowest order: abe„ — 
{iga£a + ifio-oCTfjc + -Fbe„)/F. Substitute this result into 
Eq. dHD yields abc = T~'^Vllabc - V '^ YJ^I^ ggi laotg - 

iYJ^Zi 9<7£a'5e„c, where VLq = \/j27=i ^lo- The 
Langevin noise terms are neglected in the present results, 
since under the adiabatic condition the Langevin noise 
terms have no effect on EIT quantum memory technique 



tan( 



gi92---9m-2V N 



[ET='i^%nt-Z^,9f)f' 



and 



tan (f>.i — 



Ili=i 9i^j+i,o 



[ELi{^foirsLis^^9'.)f'' 



Using the definitions above, one can transform the 
equations of motion for the probe fields and the atomic 
variables into the new field variables. With the low- 
excitation approximation and neglecting the nonlinear 
effects we find that the DSP field satisfies 
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m-2 



—9^+ ^ 't'j '^os Sj — - sin 29 cos a„i_2 9z^, 



m — 2 m — 3 



+CCOS0 ^ COS0; sin2(/ij_i( — 



1 Vj COSQfj-i. 



(10) 



where we have defined 

cosa„ = c^=^^ — „ cr = 1,2, ...,m- 3 

and — del, cos 9. It then follows that 

■Sl = nj=2^COS(/)j(-sin0i£i + C0S(/)i£2), S2 = 

nj=3^cos(/)j(-sin02(cos(?;>if ^ + sin(/)if2) + cos(/)2f3), 
and generally 

m — 3 



with 



fc fe-i 



m— 2 l—ni 



= cos(/)fc5fe+i - sin(/)fe ^ ( J]^ COS0/) sint/),,!^!^:^ 

(12) 



k-l 

sin0fc cos(/);£i. 



1=1 



Sfe = cos(j>jJk, A; = 1, 2, m — 3, (11) obtained as 

I 



On the other hand, the equation of BSP field can be 



m-2 



^^io/^^ON2V^^cos6' 



9j 



m—3 



-9f(sin6'5' — cos0<i>) + cos0(.gir2oi sin^isi — ^ g/fiio cos (/);-is;-i] 
I 



(13) 



1=2 



By comparing Eqs. (fTO| and (jT3|l with the corresponding 
DSP and BSP fields in the three-level system, one can 
see a key difference is the appearance of Sj{z,t) from 
the probe fields in our model. The adiabatic condition 
in the present case can be fulfilled only if Sj{z,t) — 
for all j. However, the input probe pulses are generally 
independent of each other so that the fields Sj need not be 
zero. To study the dynamics of the DSP field, we should 
therefore investigate first the pulse matching between all 
the probe fields needed for adiabatic condition. Bearing 



these ideas in mind, we next explore the evolution of a 
set of normal fields by introducing 

Gjj+i ^ -sin(f)j^j+i£j{z,t) + cos(l)jj+i£j+i{z,t), (14) 

where j = 1,2,..., to — 3 and tan(/)jj+i = 
gjQ,j+i^o/gj+ifljo. From the Eq. Q and together with 
the results of C7i,e^ and ai,c one can verify that the field 
Gjj+i satisfies the equation 



J 



(#iK±5|+i^cos2/3 



$12 



(15) 



r 



with 



tan^ P = 



9j+i) 



0^02 I f]2 q2 



and £j,j+i — cos (/)jj-+ifj(z, i) + sin0jj+ifj+i(z, i). 
F{£„) includes no £j or £j+i- The first term in the 
right hand side of Eq. ITSl) reveals a very large ab- 
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sorption of Gj j+i, which results in a large decay in 
the field G^j+i and then the system satisfies the pulse 
matching condition [isL 15 . 



tan 

For an explicit numerical estimation, we set some typ- 
ical values [3, gj « gj+i - 10^s~\7V « 10^ 



r 



10 s , so that the life time of field Gj 



is about < 10 s which is much smaller than the 
storage time (4|. Furthermore, by introducing the adi- 
abatic parameter r = {Y.j{'^/9jf f''^/^T where T is 
the characteristic time scale, we can calculate the lowest 
order in Eq. (fTSj) and obtain l> « 0, j+i « 0. On the 
other hand, under the condition of pulse matching, one 
can verify that Sj(z,t) cx Gj,j+i = 0. Thus equation (|10|) 
is reduced to the shape- and state-preserving case 



((9f + ccos^ cosam-2dz) ^(-z, t) = 0. 



(16) 



The formula (|T6|) is the main result of the present work. 
The group velocity of the DSP field is 



Vg = cos-" e 



E 771 — 2 



(17) 



One should bear in mind that the wave vectors kpj can 
be positive (in the +z direction) or negative (in the —z 
direction). So, by adjusting the Rabi- frequencies for ex- 
ternal pump fields properly under the adiabatic condition 
so that cosQ!„i_2 = 0, we can obtain a zero velocity for 
the DSP field. In particular, one may set No. 1 to No. 
m — 3 pump/probe pulses in the -\-z direction, while No. 
m — 2 pump/probe pulse in the — z direction (Fig. 1(b)) 



and ri„ 



-fi^g, in an experiment so that 



the group velocity = 0. In this way, we create the 
multi- frequency stationary pulses with each component: 



m — 3 



£i — COS0 Y[ cos0j4'(z,t). 



m— 3 

cos9sm(j)i-i Y[ COS (j)j'i>{z,t), 

0=1 

2, ...,m- 2. 



(18) 



These components interfere to create a sharp spatial en- 
velope. It is helpful to present a comparison of our re- 
sults with those obtained for a three- level system [^, Q : 
i) In the present system, all optical pulses can couple 
in resonance to the corresponding atomic transitions, 
thus all the applied probe fields with different frequencies 
contribute to the generation of stationary pulses. This 
means the present process is a spatial compression of ex- 
citation, which allows us to shape and intensify the lo- 
calized stationary pulses. Our technique can therefore be 
expected to enhance further nonlinear couplings and be 
applied in the most straightforward manner, e.g. to ap- 
plications in quantum nonlinear optics (l7| . Numerical 



results in Fig. 2 show how tight localization of station- 
ary pulses can be readily obtained when the multi-level 
system is applied. In contrast, for a three-level system, 
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FIG. 2: (color online) Localization of created stationary 
pulses for 5-level (red solid line) cases, where three input 
probe lights are used and the parameters are set as uje3e2 = 
uje2<ii ~ i^2/100 (a). As a comparison, blue dashed line shows 
the stationary pulses created in 3-level system (b) by employ- 
ing one standing wave of pump fields. The probe lights are 
used with the envelop exp(— z^). 



a frequency-comb is used to create a localized pulse, fil- 
tering the off-resonant input probe pulses and retaining 
only the resonant one for creation [ol]. Generally, the 
total number of probe photons created by a frequency- 
comb in a three-level system is much less than in the 
current model; ii) The present technique can be freely 
controlled. For example, based on our results, we see 
that the pulse matching in the present case is between 
all of the probe pulses with different frequencies, say, 
Ea — nj=/ tan [l > 1,(7 < TO — 2). Thus, in 

principle, one can use just one pump field to achieve the 
stationary pulse by adjusting its intensity to match those 
of the other pump fields; iii) It requires no standing waves 
in the pump fields or spatially modulated pump fields to 
create localized stationary pulses. 

Experimentally, the simplest multi-level system is an 
ensemble of four-level double A-type ^^Rb atoms. The 
schematic diagram for experimental realization is shown 
in Fig. 3. All the atoms first are trapped in state \b) 
(5^5*1/2) and only the ±z directional propagation pump 
fields {Qi and Q2) are applied to couple the transi- 
tions from |c) {5^Si/2) to |ei) {5^Pi/2{F = 1)) and {62) 
{^'^Pz/2{F = 1)) respectively. We then input the probe 
pulses (^1,2) and allow the system to achieve adiabatic 
condition. Finally, by adjusting f2i or fl2 so that gi^2a — 
g2^io, we can create the stationary pulses for probe fields 
£i{z,t) = cos 9 cos 82 {z,t) = cos sin (/)^, where 5' is 
determined by the Eq. ^ with to = 4. It can be ex- 
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FIG. 3: (color online) (a)(b) Schematic of experimental re- 
alization of stationary pulses with four-level double-A-type 
*'^Rb atoms coupled to two single-mode quantized and two 
classical pump fields that propagate in +z and — z directions, 
respectively. 



pected that when the level number m becomes larger, the 
more tightly localized stationary pulses can be created. 
According to the numerical results in Fig. 2, the effect 
becomes substantial when m > 5. 

In conclusion, we have shown the tightly localized sta- 
tionary pulses can be obtained in the general multi-level 
KIT system. We have examined the dynamics of DSPs in 
detail and found that, all the applied probe pulses with 
different frequencies contribute to the stationary pulses. 
The present process is therefore a spatial compression of 
excitation, which may be able to enhance further non- 
linear optical couplings and will have interesting applica- 
tions in quantum nonlinear optics [l7l |. In particular, our 
technique may open up a novel way towards the spatial 
compression of many probe pho tons with small losses. 
According to the results in [13|, if initially input is a 
non-classical probe pulse, e.g. a quantum superposition 
of coherent states, we may also generate entangled sta- 
tionary pulses within our model. 

We thank Dr. Y. Zhao (Heideberg University) for the 
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